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ABSTRACT 

We derive an analytic prediction for the star formation rate in environments ranging from normal 
galactic disks to starbursts and ULIRGs in terms of the observables of those systems. Our calcu- 
lation is based on three premises: (1) star formation occurs in virialized molecular clouds that are 
supersonically turbulent; (2) the density distribution within these clouds is lognormal, as expected 
for supersonic isothermal turbulence; (3) stars form in any sub-region of a cloud that is so overdense 
that its gravitational potential energy exceeds the energy in turbulent motions. We show that a the- 
ory based on this model is consistent with simulations and with the observed star formation rate in 
the Milky Way. We use our theory to derive the Kennicutt-Schmidt Law from first principles, and 
make other predictions that can be tested by future observations. We also provide an algorithm for 
estimating the star formation rate that is suitable for inclusion in numerical simulations. 

Subject headings: galaxies: ISM — hydrodynamics — ISM: clouds — ISM: kinematics and dynamics 
— stars: formation — turbulence 



1. INTRODUCTION 

The disk of the Milky Way contains ~ 10 9 M w of 
molec ular gas ( Will iams fc McKeell997HBronfman et al.l 
2000), mostly arranged in giant molecular clouds 
(GMCs) with t ypical masses of ~ 10 6 M and densities 
nn ~ 100 cm~ 3 (jSolomon et al.ll 987). Absent other sup- 
port, this gas should collapse on its free-fall time scale, 
iff ~ 4 Myr, producing new stars at a rate of roughly 
~ 250 Mq yr . However, the observed star forma- 
tion rate (SFR) in the Mi lky Way is only ~ 3 M Q yr -1 
l)McKee fc Williams! I1997D . This surprisingly low star 
forma tion rate, first pointed out bv lZuckerman &: Evansl 
remains one of the major unsolved riddles for 
theories of the interstellar medium (ISM). 

In the last 30 years, observations of star formation trac- 
ers such as Ha in other galaxies have s hown that the 
proble m is not limited to the Milky Way. Wong & Blitz 
(2002) inferred gas depletion times, defined as the ra- 
tio of the molecular surface density to the star forma- 
tion rate per unit area, of a few Gyr in resolved obser- 
vations of seven nearby galaxies. This is two orders of 
magnitude larger than the typical free-fall times of a few 
tens of Myr they inferred based on the cloud densities. 
iRownd fc Yound l|1999ft and lYoung et all (|199^ obtain 
similar gas depletion times from unresolved observations 
in many other galaxies. Nor is the problem limited to 
normal disk galaxies lik e the Milky Way. In 87 starbursts 
IGao fc Solomon! (|2004) find CO gas depletion times of 
several 0.1 — 1 Gyr, a factor of ten or less smaller than 
that in disk gal axies, and still much longer than typical 
free-fall times. iDownes fc Solomon! l)1998l) obtain rela- 
tively similar depletion times at comparable densities for 
circumnuclear starbursts in 3 nearby galaxies, and this 
range of depletion times and characteristic free-fall times 



seems typical of starbursts (jKennicuttll 1 998bl) . 

An interesting addition to this problem is that the 
star formation rate follows clear correlations. Surveys 
of many galaxies over a range of star formation rates 
and surface densities show that the star formation rate 
per unit area obeys the Kennicutt-Schmidt Law, which 
can be stated in two forms, both equally consistent with 
observations: 
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where E* is the star formation rate per unit area, E g 
is the surface density of gas, and Td yn is the dynami- 
cal (i.e. orbital) time scale of the galactic disk (Schmidt 
1959,1963; Kennicutt 1998a,1998b; Schmidt's two papers 
proposed a relationship between gas density or surface 
density and star formation rate, while Kennicutt 's deter- 
mined the exponents and coefficients of the correlations 
in equations ^ and [21 from a large galaxy sample) . Both 
forms fit the the observed sample of galaxies very well 
over a range of nearly eight orders of magnitude in star 
formation rate. 

Any successful theory of star formation must be able 
to reproduce both the lower-than-expected star forma- 
tion rate and both forms of the Kennicutt-Schmidt Law, 
and must do so using physics that is applicable in a range 
of environments from Milky Way-like disk galaxies where 
the ISM is entirely atomic and the star formation rate is 
low, to ULIRGs, where the ISM is fully molecular and the 
star formation rate is many orders of magnitude larger. 
To date, no theory is able to meet these requirements. 
Recent numerical work has been able to reproduce some 
of the observations, bu t only with cons iderable assump- 
tions and limitations. IKravtsovl 1)20031) uses the proba- 
bility distribution of densities in simulations to suggest 
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that the fraction of high density gas varies with the over- 
all density to roughly the 1.4 power, explaining one form 
of the Kennicutt-Schmidt Law. However, this observa- 
tion does not explain the other form of the Law, and it 
also says nothing about the absolute rate at which star 
formation occurs. It also fails to explain the choice of 
density cutoff that const i tutes "high density." Similarly, 
ILL Mac Low, fc Klessenl 1)200 5) show that their simula- 
tions reproduce the £ g - 4 form of the Kennicutt-Schmidt 
Law. However, their simulations depend on both an arbi- 
trarily density chosen threshold for star formation and an 
arbitrary choice of the star formation rate in gas denser 
t han the thr eshold. 

iTanl (|2OQ0fl proposes an analytic theory based on star 
formation induced by cloud-cloud collisions to explain 
the Kennicutt-Schmidt Law. In this model, the star 
formation rate is proportional to E g /rd yn because the 
inter-cloud collision time is proportional to the dynami- 
cal time, and the supply of gas available is proportional 
to the gas surface density. However, this theory also re- 
lies on an unknown efficiency of (collision-induced) star 
formation that can be roughly calibrated from obs erva- 
tions, but is not independently predicted. Similarly. iSilkl 
(1997) proposes a theory in which the star formation 
rate is set by supernova feedback. However, the theory 
depends critically on the porosity P of the interstellar 
medium to gas heated by supernovae, and it is unclear 
how P varies from the predominantly atomic, diffuse gas 
disks found in normal galaxies like the Milky Way to 
the dense, entirely molecular interstellar media found 
in starbursts. In particular, the theory predicts that, 
if P is roughly constant (as is required to obtain the ob- 
served star formation rate and the Kennicutt-Schmidt 
Law), then all galaxies should have the same ISM veloc- 
ity dispersion. This predic tion clearly fails in starbursts 
Ipownes fc Solomonlll9'9l| . 

Another broad class of theories appeals to magnetic 
fields and ambipolar diffusion. In these models, star- 
forming regions are threaded by a magnetic field strong 
enough to make them magnetically sub-critical, so the 
collapse time is set by the time required for the field to 
escape from the gas via ambipolar diffusion (see reviews 
by Shu et al. 1987 and Mouschovias 1987; for a more re- 
cent discussion, see Tassis & Mouschovias 2004 ). While 
we discuss this theory in more detail in § 17.31 we note 
that observations of magnetic field strengths in Milky 
Way GMCs, both directly via Zeeman splitting l|Crutcherl 
1999; iBourke et al.1 120011) and indirectly via statistical 
indicators l|Padoan et alJl2004[) . suggest that their mag- 
netic fields are not strong enough by themselves to pre- 
vent rapid collapse. Nothing is known of magnetic field 
strengths in other galaxies, so it is unknown if this model 
can explain the Kennicutt-Schmidt Law. 

A final class of theories, on which we shall focus, 
relies on turbulence. Observed GMCs in the Milky 
Way and in nearby galaxies have significant non-thermal 
linewidths (e.g. Fukui et al. 2001, Engargiola et al. 
2003, Rosolowsky & Blitz 2005), and this is generally 
interpreted as indicating the presence of supersonic tur- 
bulence. In a cloud supported against collapse by su- 
personic turbulence, at any given time most of the mass 
should be in structures that are insufficiently dense to 
collapse (see reviews by Mac Low & Klessen 2004 and 



Elmegreen & Scalo 2004). This conclusion is bolstered 
by simulations (e.g. Klessen, Heitsch, & Mac Low 2000; 
Li et al. 2004) that show that, under at least some cir- 
cumstances, supersonic turbulence can inhibit star for- 
mation^^ 

iPadoanl l)1995() provides an analytic theory of the star 
formation rate in a turbulent medium that depends on 
the properties of GMCs, and on the distribution of 
masses of clumps that results from turbulent fragmen- 
tation. For Milky Way GMCs it produces a value of 
the star formation rate reasonably in agreement with 
observations, but there is no way to extend this result 
to galaxie s where we c annot direct ly observe the GMCs. 
Similarly, Elm egreenl 1)20021 12003|) uses the probability 
distribution function of densities in a turbulent medium 
to estimate the mass fraction of galactic GMCs above 
a critical density of ~ 10 5 cm~ 3 , and argue that this 
can explain the low star formation rate. However, it is 
not clear why the critical density is 10 5 cm~ 3 , or how 
this value might vary from galaxy to galaxy. Nor is it 
clear how this analysis leads to the Kennicutt-Schmidt 
Law. Elmegreen argues that the law £» oc £ g - 4 can be 
explained in this picture if all galaxies have roughly the 
same scale height, but does not provide a physical reason 
why the scale height should be constant. 

Our goal in this paper is to provide a theory of the 
star formation rate that can explain both the surprisingly 
low star formation rate and two forms of the Kennicutt- 
Schmidt Law, and that can do so over a range of condi- 
tions from normal disks to ULIRGs. In other words, we 
seek to explain both the exponents and the coefficients of 
the Kennicutt-Schmidt Laws over their entire observed 
range. Our theory does not depend on an unknown effi- 
ciency or critical density for star formation. Instead, we 
proceed from three premises that are well-motivated by 
a combination of observations, simulations, and theoret- 
ical considerations. First, we assume that star formation 
occurs primarily in molecular clouds that are virialized 
and supersonically turbulent. Second, we assume that 
the probability distribution of densities is lognormal, as 
expected for supersonic isothermal turbulence. Third, 
we assume that gas collapses in regions where the local 
gravitational potential energy exceeds the local turbulent 
energy. In § [21 we develop these premises to compute the 
star formation rate in a cloud in terms of its Mach num- 
ber and virial parameter. We check this theory against 
simulations, and show that it is able to reproduce them 
well. In §|31 we apply our estimate to galaxies, and derive 
an estimate for the star formation rate as a function of 
the observable properties of galaxies. In §0J we compare 
our theoretical predictions to the observed star forma- 
tion rates in the Milky Way, and in § [3] we compare to 
a large sample of galactic-average star formation rates. 
We show that our theory provides an excellent fit to the 
data. In § we present three future observations that 
can be used to check our theory. Finally, in § [7] and § [5] 
we discuss and summarize our conclusions. 

2. THE STAR FORMATION RATE PER FREE-FALL TIME 

In this section we present our general theory of turbu- 
lent regulation of the star formation rate in dimensionless 
terms. For convenience we define the dimensionless star 
formation rate per free-fall time, SFRff , which is the frac- 
tion of an object's gaseous mass that is transformed into 
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stars in one free-fall time at the object's mean density. 

2.1. Derivation 

Both simulations and observations of turbulence in 
the interstellar medium show that the turbulent veloc- 
ity dispersion 07 computed over a volume of charac- 
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versal property of supersonic turbulence, and holds over 
a very wide range of length s cales in molecular clouds. 
lOssenkopf fc Mac Low! {2002) summarize observations 
of the Polaris Flare molecular cloud that show the 
linewidth-size relation over three orders of magnitude 
in length. Because velocity dispersions are smaller on 
smaller scales, even though the velocity dispersion may 
be supersonic over length scales comparable to the size 
of a simulation box or an entire star-forming cloud, there 
will be some smaller scale over which it is not. Vazquez- 
Semadeni, Ballesteros-Paredes, & Klessen (2003, here- 
after VBK03) show that the scale at which the turbu- 
lence transitions from supersonic to subsonic, the sonic 
length A s , is a key determinant of whether SFRff will be 
high or low. For the purposes of this paper, we adopt 
a more specific definition of the sonic length, consistent 
with that of VBK03: let 07 (x) be the one-dimensional 
velocity dispersion computed over a sphere of diameter 
I centered at position x within a turbulent medium, and 
let 

07 = (cri(x)) v (3) 

be the volume average of 07 (x) over the entire region. 
We define A s as the length I such that 07 = c s , where c s 
is the isothermal sound speed in the region. (Note that 
our A s is the same as the turbulent pressure length lp 
introduced by Wolfire et al. 2004). The linewidth-size 
relation then becomes 



07 



(4) 



While VBK03 show that the sonic length correlates 
well with SFRff, the star formation rate per free-fall 
time is a dimensionless number and the sonic length is 
a length. On dimensional grounds, there must therefore 
be another length scale that is relevant. The natural 
candidate is the Jeans length, 



At = 




(5) 



where c s is the sound speed and p is the density at a 
given point. Of course in a turbulent medium p varies 
from place to place, and we account for this effect be- 
low. Consider a "core" , a sphere of gas embedded in the 
cloud. The thermal pressure at the surface of the sphere 
is roughly pc 2 & . The largest mass such an object can have 
and remain stable against gravitational collap se is the 
Bonnor-Ebert mass llEbertlll955t IBonnorlll95^ . 
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The radius of such a sphere is roughly 

R BE « 0.37Aj. (8) 
The gravitational potential energy of the sphere is 

& 
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= -1.06 



G 3/2 p l/2 
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Here a is a geometric factor set by the sphere's mass 
distribution, and in the numerical evaluation we have 
used a = 1.2208, the value for a maximum-mass stable 
Bonnor-Ebert sphere ijMcKee fc Hollimanl [1999^1 . The 
thermal energy of the gas is 



T th = -M BE c s 2 = 1.14 \W\. 



(10) 



Using the linewidth-size relation {3J, the average turbu- 
lent kinetic energy in the sphere is 



T twh = -Mbeo- 2 (2R B e) 
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where for the numerical evaluation in the final step we 
have used p = 0.5. Thus, a Bonnor-Ebert-mass object 
has approximately equal kinetic, thermal, and potential 
energies if Aj ~ A s . If Aj ^ A s , gravity is approximately 
balanced by thermal plus turbulent pressure, and the 
object is at best marginally stable against collapse. If 
Aj 3> A s , kinetic energy greatly exceeds both potential 
and thermal energy, and the object is stable against col- 
lapse. 

Since Aj is a function of the local density, the con- 
dition Aj ^ A s for collapse translates into a minimum 
local density required for collapse. We can use this to 
compute the star formation rate, by first estimating 
what fraction of the mass is at densities higher than 
this minimum. Numerous numerical and theoretical 
studies find that the probability distribution function 
(PDF) of the density in a supersonically turbulent 
isothermal gas is lognormal, with a dispersion that 
increases with 
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Tpags^OT^azm ie^SeinadeTr _ 
Nordlund fc Padoanl Il999t lOstriker. Gammie. fc Stone 
1999HPadoan fc Nordlund! I2002|) . IPadoan fc Nordlunc 
Ij2002() find that the PDF is well-fit by the functional 
form 
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where x = p/po is the density normalized to the mean 
density in the region p^. The mean of the log of density 
is 



ln* = -f , 



and the dispersion of the PDF is approximately 
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where M. is the one-dimensional Mach number of the 
turbulent region measured on its largest scale. Let Ajo 
be the Jeans length at the mean density. The Jeans 
length at overdensity x is Aj(x) = Xjo/y/x, which we 
wish to compare to A s . We therefore define the critical 
overdensity required for collapse as 



x > x crit = <p x — — 



Ajo 

X 



(17) 



where <f> x is a numerical factor to be determined by fitting 
in § 12.21 Gas at an overdensity of a; cr i t or higher has a 
local Jeans length smaller than the sonic length, and is 
therefore unstable to collapse. The fraction of the mass 
in collapsing structures is therefore just the fraction of 
mass at overdensities of x cr ;t or greater, which is 



/ = / x*j- dx. 



(18) 



To convert / to a star formation rate, we must account 
for two factors. First, approximately 25% — 75% of the 
mass in star-forming cor es will be ejected by outflows 

We define e coro 



l|Matzner fc McKeel 1200(1 . We define e core as the frac- 
tion of the mass that reaches the collapsing core phase 
that eventually winds up in a star, and adopt a fiducial 
value of e core = 0.5. Second, we have computed the frac- 
tion of mass in collapsing structures at any given time. 
To convert this to a rate, we must divide by the charac- 
teristic time scale over which new gas becomes unstable. 
When a region collapses, it detaches from the turbulent 
flow and thereby removes pressure support from the re- 
maining, stable gas. The remaining gas will respond to 
this loss of pressure support on its gravitational collapse 
time scale, the free-fall time. We therefore estimate that 
new gas becomes gravitationally unstable over a free-fall 
timescale tg. (Alternately, we could have used a crossing 
time, which is very similar in a real GMC.) However, this 
is just a rough argument, so we let the true time scale 
be (j>ttg. We wi ll de termine <pt for purely hydrodynamic 
turbulence in § 12.21 Magnetic fields can delay collapse 
and make 4> t somewhat larger for a real cl oud than our 
fit will find ijVazauez-Semadeni et al.ll2005f) . 

With these two factors defined, the star formation rate 
per free-fall time is 



SFR ff = ■ 



xp(x) dx 



1 + erf I 
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The total star formation rate arising from a cloud of mass 
M mol is 

M* = SFR ff ^21. (21) 
tg 

We plot SFRff as a function of x cr g for <\> t — \ in Figure 
HI We can also define a "core formation rate" CFRg, 
which reflects the rate at which mass begins to collapse, 
ignoring what fraction of it will be ejected by feedback. 

T his is simply S FRff with e rn rs = 1- 

iPadoanl (|1995[) and iPadoan fc Nordlundl (|2002t I2004J) 
have previously approached the problem of estimating 
the star formation rate by considering the combination of 
the PDF of densities and the mass distribution of clumps 



-3 



_ 




1 ' ' 1 • 




■ ~- . 


- 


: \ 


X 


* 




\ 








\ \ - 

\ \ . 






\ \ 






\ \ - 


- 




\ V 






\ 






\ : 






\ : 






\ : 






\ : 






\ 






\ x ' 






\ \ : 






\ \ : 



°9 X crlt 

Fig. 1. — Star formation rate per free-fall time versus critical 
overdensity, for Mach numbers of 5 (solid line), 50 (dashed line), 
and 500 (dot-dashed line). 



created by fragmentation in a turbulent medium. Since 
we are interested only in the rate at which stars form, and 
not their mass distribution, we may neglect the clump 
mass distribution and consider only the distribution of 
densities. In so doing, we implicitly assume that all or 
most of the mass that is at densities rendering it ca- 
pable of collapse will in be in the presence of enough 
other high-density gas so that it does collapse. This as- 
sumption is bolstered by the observation that turbulence 
tends to organize mass into filaments and voids on large 
scales, so that high density gas is likely to be in the pres- 
ence of other high density gas. Moreover, as we show 
in § 12.21 this assumption produces a theory that shows 
good agreement with simulations. 

2.2. Comparison to Simulations 

To test our theory, we compare to the work of VBK03, 
who simulated a turbulent periodic box of gas and 
computed the fraction of the mass that collapsed into 
stars. The simulation setup is de scribed in detail in 
iKlessen. Heitsch. fc Mac Low! (HOQQf) , but we summarize 
it here. In simulation units, the box length is L = 2, 
the sound speed c s = 0.1, the mean density is po = 1/8, 
the Jeans length at that density is Ajo = 1/2, and the 
free-fall time is tg = 1.5. Turbulence is driven at a one- 
dimensional Mach number A4 = 2, 3.2, 6, or 10 using 
a driving field that contains power only at wavenumbers 
around k = 2, 4, or 8, where k = L/Xd and A^ is the 
driving wavelength. 

We read off the sonic length from Figure 3c of 
VBK03, noting that VBK03 define the sonic length using 
the three-dimensional velocity dispersion, while we use 
the one-dimensional velocity dispersion (J. Ballesteros- 
Paredes and E. Vazquez-Semadeni, private communica- 
tion). Given the scaling a oc l p y and assuming the turbu- 
lent velocity field is roughly isotropic, the two are related 
by A s w 3 1/(2p) A s3 . We adopt p = 0.5 through the rest 
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Run name 


A s 3 


t 


SFRff_ sim 


SFRff_th 


M2K2 


0.16 


0.62 


0.24 


0.33 


M2K4 


0.10 


0.62 


0.24 


0.18 


M2K8 


0.20 


0.62 


0.24 


0.39 


M3.2K2 


0.080 


0.44 


0.34 


0.18 


M3.2K4 


0.046 


1.58 


0.095 


0.0641 


M3.2K8 


0.031 


2.48 


0.060 


0.023 


M6K2 


0.039 


0.30 


0.50 


0.11 


M6K4 


0.023 


2.27 


0.066 


0.045 


M6K8 


0.016 


6.89 


0.022 


0.019 


M10K2 


0.018 


0.87 


0.17 


0.060 


M10K4 a 


0.013 


6.03 


0.014 


0.035 


M10K8 6 


0.0094 


4.69 


0.026 


0.018 



a Run ended with / ; 
b Run ended with / : 



0.058 of mass collapsed 
0.084 of mass collapsed 



TABLE 1 

Col. (1): Run name in VBK03. MmKfc indicates that the 
1-D Mach number is m and the run is driven at 
wavenumber k. col. (2): measured value of a s 3 . col. (3): 
Time at which 10% of the mass had collapsed, or when the 
run ended, in code units. col. (4): star formation rate in 

THE SIMULATION, DEFINED AS SFRff = //(t/tff). COL. (5): 

Theoretically estimated SFRff. 



of this section. To determine the SFRff, we read off data 
from Figure 2 of VBK03. We measure the time t at which 
a fraction / = 0.1 of the mass in the run has collapsed 
into stars. For runs where less than 10% of the mass has 
collapsed by the end, we measure / and t at the point 
where the run ends. We then compute SFRff = 1.5 f/t. 
(Using 20% instead of 10% did not substantially change 
the result.) We summarize all of this in Tabled 

We fit the VBK03 data to our theoretical estimate of 
CFRff rather than SFRff because the VBK03 simulations 
do not include any feedback. The cases with large a; cr it 
are closest to the environment in real star-forming clouds, 
so we weight by x^ rit . A Levenberg-Marquardt fit with 
this weighting gives (j> x = 1.12 and <pt = 1.91. We com- 
pare the simulation to CFRff evaluated with the best-fit 
values in Table ^ an d in Figure [5] In the Figure, the 
simulation points have error bars corresponding to a fac- 
tor of 2 uncertainty, as recommended by VBK03. As the 
plot shows, there is a large scatter, but we are able to 
reproduce the overall behavior of the VBK03 simulation 
data quite well. Note that 4>t — 1-91 implies that, for 
virialized objects, the characteristic time scale is roughly 
a single crossing time - see § 17.71 

To understand how magnetic fi elds migh t chan ge our 
results, we examine the work of iLi et all l)2004t) . who 
measure the amount of mass collapsed into cores as 
a function of time a magnetohydrodynamic periodic 
box simulation similar to those of VBK03 (identical 
box length, Jeans length, sound speed, and free-fall 
time). The initial box is magnetically supercritical, with 
M/M$ = 8.3. The simulation is driven with a three- 
dimensional Mach number of 10 (one-dimensional Mach 
number M = 5.8) at a driving wavenumber of k = 2, 
and is therefo re very similar to run M6K2 in VBK03. 
ILi et all (|2004) do not measure a sonic length, so we use 
the measured sonic length of A S 3 = 0.039 from the cor- 
responding VBK03 run. With these parameters and our 
best-fit values of 4> x and t , we find CFRff = 0.11. Read- 
ing off the time at which 10% of the mass had collapsed 
in the highest resolution run from Figure 6 of ILi et alJ 



0.10 - 




0.01 



Fig. 2. — Star formation rate per free-fall time versus Ajo/Aa, as 
measured from the VBK03 runs (error bars with diamonds) and as 
predicted by our theoretical model (asterisks). The lines show our 
theoretical predictions of SFRff versus Ajo/A s for Mach numbers 
of 2 (lowest line), 3.2 (second line), 6 (third line), and 10 (highest 
line). 



(2004) gives a measured CFRff = 0.072. The simula- 
tion result is slightly lower, but well within the factor of 
two error recommended by VBK03. While this is only 
one simulation, it provides some confidence that the in- 
clusion of magnetic fields in the supercritical regime will 
not change the star formation rate substantially. 

2.3. SFRff in Virialized Objects 

Using our theory, we can c ompute SFRff in virialized 
molecular clouds and clumps. iBertoldi fc McKeel l)1992fl 
define the virial parameter for a spherical cloud as 



5ototj 
GM ' 



(22) 



where Ctot is the one-dimensional thermal plus turbulent 
velocity dispersion over the entire cloud, R is the radius 
of the cloud, and M is the mass. Since we are concerned 
with large star- forming clouds that have a to t c s , er tot 
is approximately equal to the turbulent velocity on the 
largest scale, which we denote (J2R- Clouds with a v - n rs 1 
are in self-gravitating virial equilibrium, meaning that 
internal pressure (turbulent plus thermal) approximately 
balances gravity. Clouds with a V ir 3> 1 are non-self- 
gravitating and are confined by external pressure, while 
avir "C 1 indicates either that the cloud is supported 
against gravity by a magnetic pressure larger than either 
the turbulent or thermal pressure, or that the cloud is 
undergoing free-fall collapse. We refer to objects with 
a v j r w 1 as "virialized." 

Consider now a star-forming region that follows the 
lincwidth-size relation 



I 

2R 



(23) 



(> 
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The sonic scale is therefore 



1 000 p 



\02RJ 



l/p 



(24) 



and the Jeans length at the mean density po is 



v.to 




27TC= 



3GM 



(25) 



Thus, the critical overdensity required for collapse is 



(27) 

-^l.07M 2 , (28) 

where A4 = <J2r/c s is the Mach number of the region. 
We have used the definition of the virial parameter 12211 
in the second step, and for the numerical evaluation we 
have used our best-fit value of <p x and taken a v - a = 1.3. 
This choice is based on the evaluation o f Milky Way 
GMCs performed by McKee & Tan ( 2003). We discuss 
it in more detail in § 17.51 From this formulation, it is 
straightforward using 1(277)) to compute SFRff in a cloud 
in terms of a vlr and M. for the cloud. We have therefore 
succeeded in computing the dimensionless star formation 
rate SFRff in terms of the two basic dimensionless num- 
bers that describe a turbulent cloud: the ratio of kinetic 
to potential energy (roughly a V i r ) and the ratio of kinetic 
to thermal energy (roughly M 2 ). This relation has an 
intuitive physical interpretation. At an overdensity of 
lEcritj the thermal pressure is 

(29) 

Thus, the gas capable of collapse is simply the gas that is 
dense enough so that its thermal pressure is comparable 
to or greater than the mean turbulent pressure in the 

Cloud, Pturb 

In Figure [21 we plot the star formation rate per free- 
fall time as a function of a v - 1T and M. for p = 0.5. For 
convenience, we also fit SFRff by a power law, 



15 



Pth = RC 2 S W Pqo\ r = Pturb- 



SFRff « 0.014 



(t 



-0.68 



M 
100 



-0.32 



(30) 



Figure 0] shows the error in our power-law fit as a func- 
tion of M. and a v j r . The error is less than 5% for values 
of a V i r from ~ 0.5 — 3 and M. from ~ 10 — 1000. Since 
real star-forming clouds generally fall within this range 
(see §EJ, this power law is a reasonably good approxima- 
tion. One important thing to note about SFRff is how 
weakly SFRff varies with M. Thus, the star formation 
rate per free-fall time in a virialized cloud depends very 
weakly on the Mach number of the cloud. This is easy 
to understand intuitively. At fixed aw, increasing M. 
increases a: cr it , raising the overdensity that the gas must 
reach to collapse. At the same time, however, increas- 
ing M. increases the width of the probability distribution 
function, putting a larger fraction of the gas at high over- 
densities. These two effects nearly cancel out, which is 
why changing M. at fixed a V i r has little effect on SFRff. 

Before proceeding, we must point out one limit of our 
analysis. We have assumed that the internal structure of 



100 




Fig. 3. — Contours of star formation rate per free-fall time 
SFRff versus c? v ; r and M. The contours are labelled by value of 
log SFRff. 



100 




Fig. 4. — Contours of the error in our power-law fit for SFRff , 
defined as error = (fit — SFRff ) /SFRff . The labels on the contours 
show the value of the error. 



GMCs follows the linewidth-size relation. However, the 
OB st ar-forming clumps observed in CS bv lPlume et 
(1997) do not. They have substantially higher velocity 
dispersions than is typical for an object of their size in 
their parent GMCs, and their sizes and velocity disper- 
sions do not appear to be correlated. We interpret these 
clumps as regions of a GMC larger than a single core 
that have become gravitationally unstable and collapsed 
to higher s urface densities and pressures than the rest of 
the GMC ijMcKee fc Tanll2003|) . increasing their velocity 
dispersions. The VBK03 simulations that we have used 
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to calibrate our model do not have enough dynamic range 
to include the presence of such regions, so our estimate of 
SFRg ignores their presence. Fortunately, clumps of this 
sort constitute only a tiny fraction of the total molecu- 
lar mass of the galaxy, and are even a small fraction of 
the mass of their parent GMCs. Thus, the error we have 
made by ignoring them is negligible on the large scales 
with which we are concerned. 

3. STAR FORMATION IN GALAXIES 

In this section we will usually give surface densi- 
ties in units of M Q pc~ 2 . For convenience, we note 
that 1 M Q pc~ 2 = 2.1 x 1(T 4 g cm" 2 = 8.9 x 
10 19 Hydrogen nuclei cm -2 , and 1 M Q pc~ 2 corresponds 
to Ay = 0.045 for the local dust to gas ratio. 

3.1. The Star Formation Law for Galactic Disks 

Our formulation applies equally well to galactic disks. 
The star formation rate per unit area of a galactic disk 
is simply 



E* — 



SFRff/cMcSg 0.061 //, 



tg 



v 0.68 



(31) 



where E g is the gas surface density of the disk, /gmc is 
the fraction of it that is in molecular clouds, and tg and 
M. are the characteristic free-fall times and Mach num- 
bers in the star-forming regions of the disk. To estimate 
these quantities, we begin by considering the mean prop- 
erties of galactic disks. Note that for galaxies like the 
Milky Way, essentially all the molecular gas is in GMCs, 
so /gmc is just the molecular fraction. For starbursts, 
we also assume that all the molecular gas is collected 
into bound clouds, although this is approximate, as we 
discuss further in § 17.11 

Consider star formation in a galactic disk with a total 
surface density of E to t- The pressure at the disk mid- 
plane is then given by (cf. Elmegreen 1989 and Blitz & 
Rosolowsky 2004) 



7T 

^mp^GSgE 



g^tot 



ip^g 2 « 



^GE 2 



(32) 

where <fi mp and 4>p are constants of order unity and / g = 
Eg/Etot is the gas fraction in the galaxy. For an isother- 
mal disk consisting entirely of gas, / g = mp = <f)p = 1 
exactly. For a real galactic disk containing stars, (f>p > 1, 
because the gravity of the stars compresses the gas. We 
show in Appendix 1X1 that <fip m 3. The scale height h g of 
the gas in the disk is related to its midplane density by 

(33) 



s 2 Pg ^Gfopz 

where <r g is the gas velocity dispersion. Using these two 
expressions to solve for the midplane density gives 



Pg 



^G0pEj 
2cr2 



(34) 



To use this result, we must know cr g , which varie s from 
~ 6 km s _1 in normal disks ({Blitz fc Rosolowskvll2004|) 
to ~ 100 km s _1 in starbursts (e.g. Downes & Solomon 
1998). To estimate the velocity dispersion, we assume 
that the star-forming part of a galaxy has a flat rota- 
tion curve with velocity v mt and is marginally Toomre 



stable, so that Q ~ 1. Both assumptions are well- 
satisfied in observed galaxies ranging from normal disks 
to starbursts , and are expected on theoretical grounds 
(IQuirkl 11972 IKennicuttl fl98pt INavarro. Frenk. fc White ! 
119971 iDownes fc Solomon! 119981: iMartin fc Kennicuttl 



2001: Scliak 2002; Navarro et all 20031) . The Toomre pa- 
rameter Q is defined as l|Toomrel ll964') 



Q 



ttGE. ttGE„ 



(35) 



where k ^/2fl is the epicyclic frequency, = v TOt /r 
is the angular velocity, and r is the galactocentric ra- 
dius. We a dopt Q = 1.5 as a typical value based on the 
survey s of lMartin fc Kennicutti piOlft and lWong fc Blitd 
( 2002). However, observed Q values range from ~ 0.5 up 
to ~ 6, and spiral arms generally decrease Q. We dis- 
cuss the resulting uncertainty in the star formation rate 
in§0 

Using to eliminate er g in we ob- 

tain the mean density in a galacti c disk midplane 
l|Thompson. Quataert. & Murravll2005H . 



(36) 
(37) 



nQ 2 G 



>6.4 x 10~ 21 Q^g^o S cm " 



where il is O in units of Myr , and Qi.5 — Q/l-5 The 
corresponding free-fall time in the midplane gas is 



3^ 1/2 
32 

0.83 Qi.sOq 1 Myr 



'h , ['—) ^ /2 ^ (38) 

(39) 

Since the filling factor of molecular clouds is less than 
unity even in galaxies where the ISM is wholly molecular 
(|Rosolowskv fc Blitdl2005aD . the mean gas density in the 
star-forming clouds will be higher than this and the free- 
fall time lower. Let 4> p be the ratio of the mean molecular 
cloud density to the mean midplane density, 

^ P ^ P ~- (40) 
Pg 

With this definition, we can write the total star formation 
rate as 



^* = (J|) ^^SFRffQ-VoMcEgfi (41) 



1/2 



(42) 



■■ 0.073 M-°- 32 '^Qi.lj GMcE g ft, 



where the numerical evaluation uses our fiducial values of 
a v ir and cf>p. Noting that f2 oc T^ n , we see that our for- 
mulation already gives us something like the Kennicutt- 
Schmidt Law, equation (J5J. The Milky Way values for 
the r emaining p arameters are M ~ 25 ( Sol omon et alJ 
H981 6„ « 20 HMcKee) ll999L Q x . 5 w 1, and / G mc « 
0.25 ()Dame et al.lll987j) . which gives a numerical coef- 
ficient of 0.03 in equation 1)41(1 . with in a factor of 2 o f 
the coefficient of 0.017 determined bv IKennicuttl |[l998a) 
based on a large sample of galaxies. (Note that for the 
observational value of (f> p we are comparing the density in 
GMCs to the density in the spiral arms, which is a factor 
of ~ 4 higher than the mean ISM density - see Nakanishi 
& Sofue 2003.) Thus, our theory seems consistent with 
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the observed Kennicutt-Schmidt Law. However, our re- 
sults depend on two quantities, <p p and M., that have 
been directly observed only in the Milky Way and a few 
nearby galaxies. To completely derive a star formation 
law in terms of observables, we must compute <p p and 
M. in terms of other quantities. Fortunately, 4> p and M 
enter our prediction to the 0.5 and 0.32 powers, so we 
are relatively insensitive to errors in them. 

3.2. The Properties of Molecular Clouds 

Our goal in this section is to estimate <j) p and M. in 
terms of observables. Our strategy is to treat molecu- 
lar clouds as gravitationally bound fragments of the in- 
terstellar medium in approximate virial balance. The 
assumption of gravitational boundedness allows us the 
estimate the typical mass of GMCs, and this mass plus 
the assumption of virial balance allows us to compute 
the overdensity and velocity dispersion in GMCs. 

In the Milky Way, most molec ular gas is in clouds 
with masses of a few x 10 6 M ^ llSolomon et al J 119871 
iHever. Carpenter, fc Snellll200H). and the LM C shows a 
similar characteristic mass l|Fukui et alpOOl. The typ- 
ical m ass is somewhat low er in M33 l| Engargiola^^l.l 
12003ft and higher in M64 l|Rosolowskv fc Blitd l2005aft . 
indicating a very rough trend of increasing CMC mass 
with increasing galaxy surface density. However, all 
of these are ordinary disk galaxies, with surface den- 
sities % 100 M Q pc~ 2 . There are no observations 
that resolve individual molecular clouds in starbursts or 
ULIRGs, so we must estimate. Since GMCs appear to 
be gravitationally bound, they must have formed via 
a gravitational collapse. For this reason, their typi- 
cal mass should be roughly the Jeans mass in a galac- 
tic d i sk llKim fc Ostrikerll200H iKim. Ostriker. fc Stone] 
12001120031) . giving 



G 2 £ e 



b 

Aft 4 
•2.5 x 10 3 Qi., 



5] 3 



fin 



0! 



(43) 

(44) 
(45) 



surface pressure, so 

P ci = 0pP mp . (46) 

In an environment with a purely molecular ISM, this is 
just the ratio of the mean pressure in a gravitationally 
bound object to its edge pressure, and is ~ 2. In a pre- 
dominantly atomic ISM, cj>p is larger because molecular 
gas exists only when it is shielded by an atomic layer, and 
the weight of the bound atomic gas increases its pressure. 
We estimate (j>p « 10 — 8 /gmc, where /gmc is the molec- 
ular gas fraction of the galaxy, in Appendix [5] 

We now write down the virial theorem for a GMC, 
using a form of the theorem obtained by combining 
equation ( 24) of iMcKeel l|1999fi with equation (A7) of 
IMcKee fc Tal 120031) : 



i< 



(47) 



where S c i is the surface density of the GMC, and a v 
the standard virial parameter, 



is 



5cr 2 i-Rcl _ 

GM cl ~ GVttMciSc 



5*3 



(48) 



Equation (|47J) is quite intuitive, as it simply equates the 
GMCs internal pressure with its weight, scaled by the 
virial parameter as an indicator of how self-gravitating 
the cloud is. Together with the definition of the turbu- 
lent pressure P c \ = Pd^, (PTJl and constitute three 
equations in the unknowns p c \, a c \, and E c i. Solving for 
the molecular cloud density gives 



Pel 



375 

47T 



1/4 



p3 



1/4 



and plugging in for P c \ and M c \ gives 

. Pel /375^ 4/ ; ^ ' 




(49) 

(50) 
(51) 



where cf>p 6 = (j>p/6. The GMC velocity dispersion is 



where in the second step we have used the definition of 
Q (equation I35fl to eliminate o~ g , and £ gj 2 is £ g in units 
of 10 2 Mq pc~ 2 . In the Milky Way near the solar circle, 
where a g ss 6 km s~ x (consistent with the sound speed 
in the warm ISM - see Heiles & Troland 2003) and £ g w 
12 M Q pc -2 l|Boulares fc Coxll 1990ft . equation (jUJ gives 
Mj w 6 x 10 6 Mq. This agrees well with observed masses 
of giant atomic- molecular complexes, of which GMCs are 
the inner parts (|Elmegreenlll989(. 11994^1 . Note that the 
Toomre mass and the Jeans mass are roughly equal for 
a disk with Q 1. The Toomre mass is Mr ~ Ay£ g , 
where At ~ Qh g is the most unstable wavelength and h s 
is the gas scale height. The Jeans mass is Mj ~ h^T, g , 
so M T ~Q 2 Mj. 

Now that we have estimated the typical masses of star- 
forming clouds, we can compute their typical densities 
from knowledge of the pressures that confine them. The 
pressure at the surface of a GMC is roughly the ambient 
pressure in the midplane of a galaxy, P mp . We define </>p 
as the ratio of the mean pressure in a cloud P c \ to the 




(52) 

kms" 1 . (53) 

The numerical evaluations are for <f)p — 3 and a V ir = 1.3. 
The range of variation of </> p with /gmc is from </> p = 7.3 



for / G 



MC 



to 



2.2 for /gmc = 1. Thus, 



is 



3 — 4 times larger in normal galaxies than in starbursts. 

To convert the velocity dispersion (|52fl to a Mach num- 
ber, we must know the sound speed in the star-forming 
clouds. Observations of a galaxy can generally determine 
the temperature T in the star-forming gas, from which 
one can easily compute the sound speed c s = -^/fcsT/m, 
where m = 3.9 x 10~ 24 g is the mean particle mass, cor- 
responding to a fully molecular gas with a ratio of 10 H 
nuclei per He nucleus. However, for the purposes of nu- 
merical evaluation we can use an average sound speed. In 
the Milky Way, th e typical temperatur e in star-forming 
clouds is ~ 10 K ijSolomon et al.lll987ft . giving a sound 
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Parameter 


Value 


Qvir 


1.3 


c 


0.3 km s _1 


Ecore 


0.5 


V 


0.5 


4>p 


3.0 


<t>— 


10 - 8/gmc 


<h 


1.91 


<f>x 


1.12 


Q 


1.5 



TABLE 2 

Col. (1): Parameter. Col. (2): Adopted value. 



speed of 0.19 km s -1 . Observed starbursts have temper - 
atures in the range 29 — 46 K (Ga o fc Solomon! 12004) . 
giving sound speeds up to 0.4 km s . For the numer- 
ical evaluations in this paper we adopt an intermediate 
value of 0.3 km s^ 1 , although c s is generally directly ob- 
servable. Since the Mach number affects the star for- 
mation rate only through SFRff, and SFRff is very in- 
sensitive to Mach number, this produces relatively little 
error. We therefore estimate the typical Mach numbers 
in star-forming regions to be 



M 



■n_ j(f>pQ 2 G-E 
V2' ' 



•5.3, 



P,6 



Ql.5 1 ^g,2- 



(54) 
(55) 



(ftp 

a v ir 



6— \ 3 ^ 8 

^ Q-^FRefGMC^^) 

vir / 

(57) 



Note that while a c \ is actually the total thermal plus 
non-thermal velocity dispersion, star-forming regions are 
highly supersonic, so a cl « o- n on-thcrmai- 

3.3. The Full Star Formation Rate 

Using our calculated values for (j> p and M., the star 
formation rate per unit area of a galactic disk is 

2 19 / 8 5 3 / 8 
E * = 33/8^5/4 

.qc f j,0.34 n -1.32 1.32 v^O.68 

M Q yr" 1 kpc" 2 , 

where the numerical evaluation uses our power law fit 
for SFRff (equation |20J) and the fiducial values of all our 
other parameters, as summarized in Table |3 If one uses 
our approximation for 4>~p m terms of /gmc, this for- 
mulation of the star formation rate now depends solely 
on observables. Note that our result is diff erent than 
the stan dard scalings with E g and fl found bv lKennicuttl 
( 1998a|, and it is therefore a new prediction that can be 
tested against future observations. Also note that this 
relation should apply not just on a galaxy-by-galaxy ba- 
sis, but within an individual galaxy as well. This tOo is a 
new observational prediction. We discuss ways of testing 
these predictions in § [6] 

4. COMPARISON TO THE MILKY WAY 

We first test our theoretical prediction against the 
Milky Way. We do so in two ways to show that the our 
results are consistent. First we use the observed prop- 
erties of the molecular gas in the Milky Way plus our 
estimate of SFRff, and second we use the estimated sur- 
face densities of various MW components. 



4.1. Estimate From Observed GMC Properties 

iBronfman et al.l (|2000f) estimate that the total mass of 
GMCs inside the solar circle is M m ?] ss 10 9 M ffl . The 
mass distribution of the clouds is iWilliams fc McKeel 
fl997h 



dN 



d In M c 



0, 

1QM 



cl.fi 



M clfi > 6 
M c , 6 < 6 



(58) 



where M r i is the cloud mass and M c \ t e — M c i/(10 6 M Q ). 
iSolomon et~aTl l|1987fl catalog 273 galactic GMCs ob- 
served in CO They find that the average column density 
of GMCs is A^h ~ 1.5 x 10 22 cm~ 2 independent of mass, 
where the subscript H indicates that we are re ferring to 
the number of hydrogen nuclei. iMcKed l)1999fl uses this 
result to estimate that the free-fall time in a GMC is 



iff = 4.7 



M r - 



10 6 Mr. 



1/4 



Myr. 



(59) 



Combining the linewidth-si z e and mass- radius re l ations 
inferred by ISolomon etaf] l)1987|) and IMcKed (|1999f) 

gives a Mach number-radius relation 



Md = 25 



\10 6 M G 



0.25 



(60) 



for a GMC temperature of 10 K. From these relations, 
it is straightforward to estimate the total star formation 
rate by integrating the star formation rate over the GMC 
mass distribution, 



6xl0 6 M 



SFRff (Md) dN 



'10 4 M Q 

w5.3 M yr -1 . 
have imposed a lower 



tff(M cl ) d\nM c 



m nn,i= / ' , ;• dM c i(6i) 

(62) 

We have imposed a lower cutoff of 10 4 Mq be- 
cause a vir ^> 1 for GMCs with smaller masses 
ijHever. Carpenter. &: Snellll200Hh which greatly reduces 
their star formation rate. The observed star forma- 
tion rate in the Milky Way is M* w 3 M yr -1 
ijMcKee fc Williams1ll997j) . so our estimate agrees with 
observations to a factor of 1.8, a reasonable fit. 

An important subtlety of this analysis is that we must 
impose a lower cutoff when integrating the star forma- 
tion rate over the GMC mass distribution because small 
clouds, if they are virialized, contribute significantly to 
the star formation rate. The integrand in l|61|l scales 
as roughly MT : one gets an exponent of —0.6 from 
the logarithmic mass spectrum dN /dlnM c i, —0.25 from 
the free-fall time, and ~ —0.08 from the dependence of 
SFRff on the Mach number, and hence on M c \. Thus, 
each decade range in the mass of virialized clouds con- 
tributes almost equally to the star formation rate. The 
contribution to the total star formation rate from small 
clouds is small not because the clouds contain a small 
amount of mass, but because small clouds are not virial- 
ized. 

4.2. Estimate From Surface Densities 

We can also compute the Milky Way star formation 
rate using surface densities, the rotation curve, and the 
velocity dispersion. The vast majority of star formation 
in the Milky Way occ urs in a ring from 3 to 11 kpc in 
galactocentric radius (jMcKee fc Williams! 119971) within 
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which th e molecular and ato mic gas surface densities are 
roughly (Wolfire ct al. 2003) 



^mol 



and 



6.3 exp ^— 
4.1 exp (- 



(rfc-4.8.5)' 
2-2. 25 a 

r k -6.97 \ 
2.89 I 



M pc" 
M pc- 2 , 



3 < r k < 6.97 

r k > 6.97 
(63) 



^HI 



(2.0r fc - 0.8) Mq pc" 
7 Mq pc" 2 
[-1.57 + 8.57(r fe /8.5)] 



r k < 4 
4 < r k < 8.5 
M pc- 2 , r k > 8.5 
(64) 

where r k is th e galactocentric radiu s in kpc, and we have 
multiplied the Wol fire et al.l l)2003|) values for the surface 
density of hydrogen by 1.4 to get the total surface density 
including both H and He. From these surface densities 
we can directly compute £ g , /gmc, and 4> P - The galactic 
rotati on speed is v rr ,t ~ 220 km s _1 , and is flat over the 
ring (Binnev & Mcrrificld 1998), so 



n 



0-22 „ i 
Myr . 

r k 



(65) 



. We estimate M. as a function of radius from S g and 



The temperature in the molecular gas is ~ 10 K 
ijSolomon et al.lll98^) . giving a sound speed c s « 0.2 km 
s 

f2 using (|5l} 

The final step is to estimate Q, which we do in two 
different ways. First, we compute Q from E g assuming 
that the gas velocity dispersion is <7 g = 6 km s _1 inde- 
pendent of radius. Th is is consistent with observations 



(|Kennicuttl 119891: iHeiles fc Trolandl 12003). although the 
observations are quite uncertain because it is difficult to 
determine the velocity dispersion as a function of radius 
within the galaxy. Second, we compute Q from the gas 
scale height, which can be directly measured in the Milky 
Way. Equation (|33|l allows us to compute p g from E g and 
h s , and equation i|36|) gives Q in terms of p g and Q. The 
scale heights of the at omic and molecular gas within the 
star-forming ring are ijWolfire et al.l l2003) 



65 pc, 

65 exp[(r fc - 8.5)/6.7] pc, 



r k < 8.5 
r k > 8.5 



and 



(66) 



33 pc, 
33 exp [(r k 



8.5)/6.7] pc, 



r k < 
r k > 



(67) 

Note th at we have converte d the half-density heights 
given bv I Wolfire et~aTI (Hki) to scale heights by assum- 



ing an isothermal density profile p oc sech [z/(2h g )]. We 
determine a Q by computing the midplane density of 
atomic and molecular gas, and then solving (|36|l for Q us- 
ing the surface density- weighted average of the two mid- 
plane densities. The result agrees to within 20% with the 
value Q as a function of radius we derive using the first 
method. We plot the azimuthally-averaged Q versus ra- 
dius for the Milky Way in Figure|5| However, most Milky 
Way s tar formation occurs in the spiral arms. Balbus 
( 1988) shows that the local Q value in a spiral arm is 
related to the azimuthally averaged Q by 



Qa 



-1/2 



(68) 




(kpc) 



• 2 



Fig. 5. — Predicted value of Q versus radius, estimatated using 
azimuthal averages and scale heights (dot-dashed line), using az- 
imuthal averages and <r g = 6 km s — 1 (dashed line), and corrected 
for spiral structure (solid line). 



The Mi lky Way's spiral arms ar e overdense by factors 
of ~ 4 (INakanishi fc Sofuei r2003). so we reduce our es- 
timated value of Q by a factor of 2 to account for this 
effect. We also show the corrected Q in Figure [SJ 

Integrating over the star-forming ring, we find a pre- 
dicted star formation rate 

/>11 kpc 

7.V ^ / nr f .0.34 n -1.32 1.32 V 0.68 

JW*-prcd~ / 9.5 JGMC <Pj; q Ql.5 "o ^g,2 

J 3 kpc 

2tt RdR Mq yr- 1 kpc~ 2 (69) 
«4.5 Moyr- 1 . (70) 

This agrees with the obser ved star formation rate of 3 M Q 
yr _1 in the Milky Way l|McKee fc Williams! Il997|) and 
with our estimate based on observed GMC properties to 
better than a factor of 2. If we omit the correction for 
spiral arms, we get a star formation rate of 2.1 Mq yr -1 , 
still in good agreement, so the spiral arm correction is 
not critical. 

Note that ^ gives a prediction not just for the to- 
tal star formation rate in the galaxy, but also for the 
radial distribution of star formation. We show this in 
Figure El For compariso n, we also show the model of 
iMcKee fc Williamsl l)1997fl (scaled to have the same in- 
tegrated star formation rate as ours), which is generally 
consistent with observational data on the radial distribu- 
tion of star formation outside 4 kpc. Our model is similar 
to the McKee & Williams model in this range, but differs 
substantially inside 4 kpc because McKee & Williams use 
a simple exponential distribution with a cutoff for the ra- 
dial variation of the molecular gas surface density, while 
we use a more accurate distribution that better reflects 
the decline in the molecular gas surface density in the in- 
ner galaxy. We find that the characteristic radius of star 
formation in the Milky Way, defined as the radius within 
which half the star formation occurs, is i? c har ~ 7.1 kpc. 
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Fig. 6. — Predicted variation in the star formation rate per unit 
area, S*, with galactocentric radius r. The solid line is our mode l, 
and the dashed line is the model of McKcc & Willi ams! (119971) , 
scaled to have the same integrated star formation rate that we 
predict. 



Taking the outer radius of the star forming disk to be 
11 kpc, this gives i? c har = 1.3i?i/ 2 . This suggests that 
the common observational practice of measuring quanti- 
ties such as ang ular velocities at h alf the outer radius of 
star formation ()KennicuttJll998bj ) should be reasonably 
accurate. 

5. COMPARISON TO GALACTIC-AVERAGE STAR 
FORMATION RATES 

5.1. Statistical Comparison 

For a second test we compare our prediction against 
a sample of 95 galaxies, ta ken from the norm al galaxies 
and starbur sts compiled by [K ennicutt ( 1998a) plus star- 
bursts from lDownes fc Solomon! l)1998|) . For the Kenni- 
cutt galaxies, we use the measured values of £ mr >i, Sp., 
and TH v r i as reported in Tables 1 and 2 of IKennicutti 
( 1998a) to compute a theoretical star formation rate from 
lfSl)|) . We follow Kennicutt in taking ft — 47r/rd yn to be 
the typical value of f2 in the star-forming region, and we 
exclude galaxies for which there is no measured value of 
Td yn . For starbursts where there is no measured value of 
E g (only S m pi) we assume /gmc = 1- 

For the IDownes k, Solomon! l|1998|) sources, we 
use a compilation of support i ng in formation from 
iThompson. Quataert. k, Murray! l)2005|) . As with the 
Kennicutt starbursts, we take /gmc = 1 for all these 
points. We derive £ g and Td yn from the gas mass, half- 
pow er radii, and rotation curv es from Tables 4, 5, and 
9 of IDownes fc Solomon! l)1998j) . and we derive star for- 
mation rates from the FIR luminos ities taken from the 
texts of IDownes fc Solomon! llT998l) (TRAS00057+4021 
IRAS02483+4302, VII Zw 31). iGenzel et alJ J200ll 
(IRAS23365+3604, IRAS17208_^014) i Ifleckmanet 
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(2000) (1RAS10565+2448L and lSoifer et al.H200ffl CMrk 
231). We compute the star formation rates from the 



FIR lum inosities using the conversion factor o flKennicuttJ 
(1998b). The data set includes mul tiple points for Arp 
193, M rk 273, and Arp 220 because IDownes fc Solomon! 
(1998) break the sources up into a more diffuse compo- 
nent and one or two "extreme" starburst nuclei. For 
these objects we include both the diffuse component 
and the nucleus or nuclei. Data for the surface densi- 
ties, dynamical times, and luminosities for the diffuse 
co mponents come from Table s 4, 5, and 9 and the text 
of|D ownes fc Solomon! l)1998[K while data for the nuclei 
come from Table 12. 
To compare to this sample, we compute 



,2 



Y 



A da 



A fi , 



-data 

) - log(XU_ thoory )] 

(71) 



for our model, and, as a normalization, for the lKennicuttJ 
(1998a) empirically determined best fit. The number of 
fit parameters JVjat is unity for the Kennicutt best fit and 
zero f or our model. We find x 2 = 0-40 for the best- 
fit of Kc nnicuttl l)1998alk while our theoretical model 
gives x 2 — 0.55. Note that these are not traditional 
X 2 goodness-of-fit statistics, since we are using the loga- 
rithm of the data, we have no error bars for the measure- 
ments, and the dominant errors (arising from extinction, 
an imperfectly known IMF, and similar astrophysical un- 
certainties - see Kennicutt 1998b) are systematic and 
therefore highly non-Gaussian. Instead, the meaning of 
this statistic is that 10 x is the RMS factor by which the 
model errs in estimating the star formation rate. Thus, 
our results corresponding to RMS errors of a factor of 
5.6 for our model and a factor of 4.3 for the Kennicutt 
fit. Given the factor of several systematic uncertainties 
in the measured star formation rates, these values are 
essentially identical. 

5.2. The Kennicutt- Schmidt Law 

The Kennicutt-Schmidt Law correlates the star forma- 
tion rate with either E g f2 or £ g , while our theory makes a 
prediction based on E g , f2, and /gmc- From an intuitive 
physical standpoint, one would be surprised if the star 
formation rate did not depend on all three of our param- 
eters to at least some degree. Thus, the two forms of the 
Kennicutt-Schmidt Law represent two ways of project- 
ing a four-dimensional space (consisting of E g , fl, /gmc, 
and £*) onto two dimensions. To compare our theory 
directly to these laws, as opposed to the underlying data 
as we did in § 15.11 we must make some additional approx- 
imations. We stress that we make these approximations 
only for the purposes of the projection, and that the right 
way to test our theory is to use the measured values of 
E g , fl, and /gmc, as we did in § 15.11 We make them 
because they allow us to show, in a relatively intuitive 
manner, why projecting the four-dimensional data down 
to two-dimensions still allows such a good fit to the ob- 
servations. 

Since neither version of the Kennicutt-Schmidt Law in- 
volves /gmc , we must estimate it in terms of or Q, . 
IWong niitl l|20fil and IRosolowskv Rr. Rlitd l|2005hfl 
find that the ratio of molecular to atomic gas follows 
the approximate relation 



J mol 
^HI 



2.5 x 10 4 cm" 3 K 



1.0 



(72) 
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with about half a dex of scatter. Since P mp is just a 
function of E g in our model, equation (|72J) gives 

/gmc- (1 + 0.025 E-j)" 1 , (73) 

for our fiducial <j)p = 3. Note that most of the dynamic 
range of the Kennicutt-Schmidt Law lies above E g = 

100 M Q pc~ 2 , for which /gmc ~ 0.98, where have made 
the approximation that most molecular gas is in GMCs. 
Thus, /qmc is almost constant over most of the range 
of the the Kennicutt-Schmidt Law, which is part of the 
reason that a projection of the data that neglects /gmc 
makes little difference. 

With /gmc approximated in terms of E g , the remain- 
ing step is to project from the three-dimensional space of 
£*, Eg, and Q onto a two-dimensional space of E* and 
just S g or E g f2. To do this, we make use of the fact that 
Eg and Q for galaxies appear to be correlated, as shown 
in Figure [7J The correlation is fit reasonably well by the 
rule 



fl Q = 0.058 E 



0.49 



(74) 



as the Figure shows. We can use this rule to estimate E g 
and f2 independently from any combination of E g and f2, 
allowing us to project our theory into the same lower- 
dimensional space as the Kennicutt-Schmidt Law. This 
correlation is the other half of the reason that project- 
ing the data into two dimensions works well: E g and Q 
are not really independent, at least in the available data 
set. Because they are correlated, projecting the data onto 
any appropriately chosen combination of them will work, 
which is why the Eg 4 and E g f2 forms of the Kennicutt- 
Schmidt Law work equally well, as does our prediction, 
which is approximately E* cx Eg 68 i7 132 . Even though 
the parameter space is four-dimensional, most of the data 
points lie near a line within it, which makes distinguish- 
ing different models quite difficult. We discuss how to 
break this degeneracy in § El Also note, however, that 
while (|74JI holds between galaxies, it is unknown if it 
holds within galaxies. For this reason, the projection we 
derive to compare to the Kennicutt-Schmidt Law may 
apply only to averages over many galaxies, not within 
individual galaxies. 

Using equations 1(751) and l(7l^l in equation JSHJ), our 
theoretical prediction for the star formation rate in terms 
of £!E„ is 



£* « 3.2 <pg 



0.34 r) -1 - 32 
P,6 "*1.S 



/gmc (^oE g , 2 ) ' 89 M yr 1 kpc 



„-l 1, -2 



where 



/gmc « 1 + 5.5 x 1(T 3 (fi E g , 2 ) 



-1.34 



(75) 
(76) 



and <j>p sa 10 — 8 /gmc- The observed Kennicu tt-Schmidt 
Law w ith this choice of dependent variable is ( Ke nnicuttl 
Il998aft 

E* = 0.017 ftEg. (77) 

We plot this and our theoretical prediction in Figure H3 
As the plot shows, our theoretical prediction, when we 
take into account the way that /gmc, E g , and are 
related, essentially reproduces the first form Kennicutt- 
Schmidt Law. If we instead choose E g to be our in- 
dependent variable, following the second form of the 
Kennicutt-Schmidt Law, our theoretical prediction is 

E* = 0.19 4 3 6 4 Q:i- 32 /cMcE g f. (78) 




2 3. 
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Fig. 7. — S g versus Q for observed galaxies. The data points 
are: normal disks (crosses), circumnuclear starbursts (asterisks), 
ULIRGs (diamonds), Arp 193 (triangles), Markarian 273 (empty 
squares) , and Arp 220 (filled sq uares ) . For a description of how we 
derived these data points, see S I5.1I The line shows a linear fit to 
the data. 



where /gmc is given by equation i|73|) and <f>p is approx- 
imated in terms o f /gmc as for the previous case. The 
observed law is is 



E* = (2.5 ± 0.7) x 10~ 4 
M Q yr -1 kpc~ 



1 M pc~ 



1.4±0.15 



(79) 

w 0.16 E g ; 4 M yr" 1 kpc" 2 (80) 

We plot this and our theoretical prediction in Figure El 
and find that, again, our fit is reasonably good. The only 
exception is at values of E g ^ 10 M pc -2 . The error 
there arises from the fact that almost all the galaxies with 
Eg so small lie well above the fi versus E g correlation 
we have used to project our theory (as shown in Figure 
01, so the values of VL we are using are systematically 
smaller than those of the real galaxies in that region of 
parameter space. Since our star formation rate depends 
on fi 1 - 32 , this underestimation of £1 causes the theory 
to underpredict the star formation rate. If one uses the 
measured values of Q rather than the linear fit, the error 
at small E g is no larger than it is elsewhere. 

6. FUTURE OBSERVATIONAL TESTS 

Our theory makes three observational predictions that 
should be directly testable in the next few years. First, 
we can test our theory on nearby galaxies where molec- 
ular clouds are directly observable. In § 14.11 we compute 
the star formation rate in the Milky Way by integrat- 
ing over the observed distribution of Milky Way GMCs. 
While we have some information about larger GMCs in 
nearby galaxies, small GMCs make a non- negligible con- 
tribution to the star formation rate there just as they do 
in the Milky Way. To reliably compute the star forma- 
tion rate in another galaxy, we must therefore identify the 
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Fig. 8. — Predicted star formation rate versus E g /T(j yn (solid 
line). We also plot the Kennicutt (1998b) best fit (dash e d line) , 
and observed points for normal galaxies from IKe nnicuttj 11998a) 
(crosses), circumnuclear starbursts from IKennicuttl 11998a!) (as- 
terisks), ULIRGs (diamonds), Arp 193 (triangles), Markarian 273 
(empty squares), Arp 220 (filled sq uares ). For a description of how 
we derived these data points, see S I5.1I 




log E g (M e pc- ! ) 

Fig. 9. — Pr edicted star forma tion rate versus S g (solid line). 
We also plot the lKennicuttl I1998bl) best fit (dashed line). The data 
points are observed galaxies: normal disks (crosses), circumnuclear 
starbursts (asterisks), ULIRGs (diamonds), Arp 193 (triangles), 
Markarian 273 (empty squares), and Arp 220 (filled sq uares ). For 
a description of how we derived these data points, see S I5.1I 



lower mass cutoff below which molecular clouds become 
non-virialized. This cutoff has not yet been observed in 
any galaxy but the Milky Way, but such an observation 
is a straightforward extension of existing data sets to 
higher sensitivities and angular resolutions. It should be 



within the capabilities of the SMA, CARMA, or ALMA. 
Once one has determined the full cloud mass distribution 
for another galaxy down to the non-virial cutoff, one can 
compute the star formation rate in another galaxy by us- 
ing H21J) to compute the star formation rate in each cloud, 
just as we have done for the Milky Way. Since this type 
of test depends only on our calculation of SFRg , and not 
on any of our calculations of GMC properties in external 
galaxies, this method allows our theory of SFRg to be 
tested independently of the rest of our model. 

Observations that resolve the star formation rate in 
a galaxy in annular rings, but cannot resolve individual 
GMCs, provide a second possible test of our theory. With 
sufficiently good data, one could use equation l|57l) to pre- 
dict the star formation rate versus radius within a galaxy, 
just as we have done for the Milky Way in § 14.21 This 
could then be compared to resolved obser vations of star 
forma tion versus radius, similar to those of lWong fc Blitd 
(2002). The primary observational challenge in this com- 
parison is that, to compare to a single galaxy where one 
cannot assume that parameters such as Q have their av- 
erage values, one must measure all the information that 
we measured for the Milky Way. In particular, one must 
know Eg, /gmCj ^, and Q as a function of radius. The 
first three are relatively straightforward, but measuring 
Q requires that one be able to measure either the ve- 
locity dispersion or the gas scale height. Neither quan- 
tity is easy to determine observationally, but without it 
the theoretical predictions will be uncertain by a factor 
of several. We suggest two possible ways to make this 
measurement. First, one could perform resolved obser- 
vations of a starburst galaxy, where <r g is large enough 
to be comparable to the galactic rotation velocity and is 
therefore easier to measure. Second, one could measure 
oy in a norma l disk that is face-on, and then use the 
Tullv & Fisher] l)1977l) relation to obtain a rotation ve- 
locity Since there is some scatter in the Tully-Fisher re- 
lation, this procedure would likely need to be performed 
over several galaxies to minimize the errors arising from 
the uncertainty in the rotation curve. 

A third possible test involves usin g a sample of galax - 
ies similar to but larger than that in IKennicuttj (^998a). 
We found in § 15.21 that, because E g and fl are them- 
selves correlated, E* will correlate equally well with an 
infinite number of combinations of E<r and Q. Our 

theory predicts that the true scaling should be E* oc 
/GMcEg 68 ^ 1 ' 32 , but the current data cannot distinguish 
this combination of E g and Q from any other. However, 
there is no reason that a future, larger sample of galax- 
ies could not. In order to break the degeneracy between 
combinations of E g and f2, a future sample must contain 
a large number of galaxies, or annuli within galaxies, with 
fixed Eg and varying Q, or fixed f2 and varying E g . With 
such a sample, one could compute predicted star forma- 
tion rates using (|57|l and repeat our analysis in S 15. H and 
determine whether Eg 68 ^ 1 - 32 is a better fit. However, 
there is likely to be considerable scatter arising from the 
stochastic nature of the cloud and star formation process. 
This test will therefore require a considerably larger sam- 
ple of galaxies than are currently available. 

Finally, note that one cannot easily test our theory by 
looking at individual GMCs. Simulations of turbulence- 
regulated star formation show significant fluctuations in 
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the star formation rate versus time, and we expect that 
real GMCs will also have large fluctuations. Thus, our 
theory is valid only as an average over an ensemble of 
GMCs. Furthermore, observations of a single GMC run 
into a problem with GMC ages. Tracers of star formation 
such as FIR and radio continuum luminosities measure 
the mass of stars formed over some period in the past. 
The amount of time depends on the tracer, but even trac- 
ers sensitive only to the very youngest stellar populations 
integrate the star formation over several Myr. We do not 
know the GMC lifetime well, and even in our model of 
virialized GMCs we cannot rule out the possibility that 
it is only ~ 10 Myr, a few GMC crossing times. Thus, 
one cannot be confident when observing a single GMC 
that it has been forming stars for long enough to have 
reached a steady-state luminosity in the tracer that one 
is using. This makes observations difficult to interpret, 
because one cannot break the degeneracy between the 
star formation rate and the age of the cloud. 

7. DISCUSSION 

7.1. Estimate of Uncertainties 

We begin to estimate our uncertainties by considering 
how much our estimates of the star formation rate could 
be off by considering a "worst-case scenario" for our un- 
known parameters, a V i r , 4>p, 4>7i Qi an d £core- Our fidu- 
cial value for 1.3, and a plausible range based on 
the observations is 1 — 2. As discussed in Appendices lAl 
and[Bl the plausible ranges in <pp and (j>p- are <f>p = 1 — 6 
and <pp = 2 — 10. We have also used a fiducial value of 
Q = 1.5. Simulations of purely gaseous magnetized disks 
show that collapse in a disk can set in at Q in the range 
0.9- 16 llKim fc OstrikerlEjOlt iKim. Ostriker. fc Stond 
120021 120031) . Analytic work sh ows that stars allow 
smaller values of Q to be stable ijJog fc Solomon! Il984t 
Rafikovl |2001[) . Observatio nally, most galaxies fall in 
the range Q = 75 - 3 l|Martin fc Kennicuttl I200H 
IWong fc Blitzi200"2T) . with outliers going as far as Q = 0.5 
to Q = 6. We adopt Q — 0.75 — 3 as our plausi- 
ble range of variation for most galaxies. Finally, we 
have taken e coro = 0.5, but the plausible range for the 
mass fraction ejected by feedback is e corc = 0.25 — 0.75 
l|Matzner fc McKedl200q) . 

If we consider all of these parameters simultaneously 
assuming their extreme values, for a given galaxy we can 
reduce our predicted star formation rate by as much as 
a factor of 10, and increase it by as much as a factor of 
6, relative to our fiducial case given by the parameters in 
Table |21 A more realistic estimate of the error is proba- 
bly a factor of ~ 3, because there is no reason our errors 
should add up systematically in this fashion. Indeed, the 
maximum errors are possible only for combinations of pa- 
rameters that can be ruled out on observational grounds 
other than the Kennicutt-Schmidt Law. For example, a 
reduction of the star formation rate by a factor of 10 is 
possible only if <j>p — 2, 4>p — 6, and a vlY = 2, giving 
4> p = 1.3. This corresponds to a galaxy where molecular 
clouds are only overdense relative to the mean in the ISM 
by 30%. No known galaxy, including ones where the ISM 
is entirely molecular, has clouds with such a small over- 
density compared to the rest of the ISM. Indeed, such a 
galaxy would effectively have no clouds at all, just a con- 
tinuous intercloud medium. Similarly, an increase in the 
star formation rate by a factor of 6 occurs for (f>p = 10, 



4>p = 1, and a v ir = 1- Plugging in Milky Way values 
of ft and Eg with these parameters gives M ~ 5, much 
smaller than the observed velocity dispersion in GMCs 
in the Milky Way or in any other galaxy. 

We can also identify a number of uncertainties 
not associated with any specific parameters, but in- 
stead with conceptual assumptions that we have made. 
First, observations have confirmed that, in at least 
some clouds, the star format ion rate is lower i n the 
outer than the inner parts llLi Evans, fc Ladal 119971 
Uohnstone. Di Francesco, fe Kirkl 120041) perhaps due to 
increased ionization there l|McKeel l1989V The periodic 
box simulations we have used to calibrate SFRg do not 
include any effects arising from the finite size of real 
GMCs, and this may produce an error. A second effect 
is that we have assumed that all the gas in starbursts is 
in bound structures capable of forming stars. However, 
observed galactic nuclei and starbursts that are molec- 
ular throughout consist of a coll ection of cloud s with 
a molecular intercloud m edium ijSolomon et alJ 119971 
iRosolowskv fc Blitzl l2005aj) . Our assumption that all 
the gas is in bound structures may therefore cause us 
to syste matically overestimate the st ar formation rate. 
However, IRosolowskv fc Blitzl i|2005a|) find that in M64 
the clouds account for ~ 75% of the mass, and simula- 
tions such as those VBK03 of show that 50% of the 
mass does collapse in unstable environments, so the er- 
ror is probably small. Third we have neglected magnetic 
fields. We argue in § 17.31 that star- forming clouds are 
likely magnetically supercritical, and thus cannot be held 
up against collapse by magnetic fields, and we present 
preliminary evidence in § 12.21 that a magnetic field in a 
fairly supercritical cloud does not substantially inhibit 
star formation. However, it is poss ible that a ma gnetic 
field stronger than the one used in iLi et all (|2004fl . yet 
still not strong enough to make the cloud subcritical, 
could inhibit the formation of cores by preventing gas 
from flowing across field lines to accrete onto them. Wc 
find this unlikely, however, since in a supercritical cloud 
the Alfven Mach number is likely to be unity or greater. 
Fourth, we have ignored the possible effects of star forma- 
tion in objects like the clumps observed bv lPlume et alJ 
( 1997) that do not lie on the linewidth-size relation, and 
that numerical simulations thus far lack the resolution to 
model. Since these objects are over-pressured and over- 
dense compared to typical galactic star- forming clouds, 
they have shorter free-fall times and form stars faster. 
By neglecting them, we probably systematically under- 
estimate the star formation rate. The extent of the un- 
derestimate is somewhat uncertain, since simulations to 
date have not modeled this effect, and we do not know 
how much exactly mass is in these clumps in the galaxy. 

7.2. Application to Simulations 

Our theory of turbulence-regulated star formation is 
readily applicable to simulations on cosmological or 
galactic scales that do not have enough resolution to 
model molecular cloud formation or star formation di- 
rectly. This is particularly true because, while we can 
integrate over an entire galactic disk to compute average 
star formation rates, we also predict the star formation 
rate in terms of local properties of the gas. 

In a simulation, one usually wants to implement star 
formation as a sub-grid model. This requires a recipe for 
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determining at what rate the mass in a given cell or par- 
ticle is transformed into stars. Equation (|21J) gives the 
star formation rate in terms of the local free-fall time iff, 
molecular mass M mo i, and the star formation rate per 
free-fall time SFRff, which is a function of a v ir and M., 
the local virial parameter and Mach number. Since in a 
simulation the density of every cell is generally known, 
it is simple to compute tg. Since the gas mass but not 
the molecular mass of every cell is known, one must de- 
termine /gmc to find M mo \. To do this, one may either 
assume that sufficiently dense cells are entirely molecu- 
lar, or more directly use the observed relation between 
pressure and /gmc, equation Ij72(l . (One should be wary 
of applying this rule to galaxies with metallicities too 
different from that of the Milky Way, though, since the 
correlation almost certainly has some metallicity depen- 
dence.) 

Finally, to compute SFRff , one needs to know the virial 
parameter and Mach number within a cell. The easiest 
way to estimate M is to compute the velocity dispersion 
over a small region around the cell, and extrapolate down 
to the size of the cell using the linewidth-size relation 
cr cx I - 5 . This plus the temperature of the cell yields M. 
within the cell. While this procedure is somewhat uncer- 
tain because it requires extrapolation to scales below the 
grid size, M. has only a weak effect on SFRff. One can 
compute a V h- by using M. to estimate the kinetic energy 
in the cell, and comparing to the estimated gravitational 
self-energy of the cell. This process for estimating a v i r is 
similar to the process of estimating whether a region is 
bound us ed in the sink particle c r eation procedures out- 
lined by iBate. Bonnell. fe Price! lll995l) for L a grangian 
codes and iKrumholz. McKee. fc Kleinl l|2004D for Eule- 
rian codes. From a v i r and M, one can compute SFRff, 
and from that M*. 

This procedure provides a simple estimate for the rate 
at which a cell turns its mass into stars that is based on 
a physical model rather than an arbitrary density cutoff 
and efficiency for star formation, which are commonly 
used in simulations now. One caveat on our approach, 
however, is that it does not apply to primordial star for- 
mation, where the primary limit on star formation is the 
ability of the gas to cool, rather than turbulent support. 

7.3. Magnetic Fields 

Our theory of star formation regulated by supersonic 
turbulence is only valid if star formation occurs primarily 
in regions that are magnetically supercritical. If molec- 
ular clouds are magnetically subcritical, then magnetic 
fields can prevent collapse, and the time required for the 
flow to "replace" collapsing cores is the ambipolar dif- 
fusion time rather than the free-fall time. This effect 
could inhibit flow in unbound regions of GMCs even if 
the clouds overall are su percritical . How ever, our com- 
parison with the work of lLi et all (|2004Tl gives prelimi- 
nary evidence that this effect is small. 

Both observations and general theoretical considera- 
tions support the idea that molecula r clou ds are super- 
critica l. Theoretically, iMcKeel l|1989l) and iMcKee et all 
( 1993) point out that GMCs cannot be bound, turbulent, 
and magnetically subcritical. Turbulence and magnetic 
fields together can support a larger cloud mass than mag- 
netic fields or turbulence alone. If the turbulent energy 
is comparable to the magnetic energy, as both observa- 



tions and general expectations of equipartition suggest, 
then the critical mass arising from both sources must be 
M cr it ~ 2M$. If the cloud is magnetically subcritical, 
then M < M$, so M < M crit /2. However, for a cloud to 
be bound it must be near its critical mass. A cloud that 
is only half its critical mass must be unbound, and would 
certainly not be centrally concentrated. Since observa- 
tions indicate that molecular clouds are both bound and 
centrally concentrated (see § I7.5|) . it follows that they 
must be magnetically supercritical, with M ss 2M$. 

Zeeman splitting observations of magnetic field 
strengths in Milky Way GMCs supp ort this view. 
iCrutcherl lH99l a n dlBourke et alJ ll200l find tha t M w 
2M... lOalli et all l|1999ft and lAllen fc ShiJ 1)200(1) point 
out that this conclusion depends on the assumed cloud 
geometry along the line of sight, and that Crutcher's data 
are consistent with M w M$ if clouds are highly flat- 
tened. However, this model is feasible only if true mag- 
netic field values in regions with no detec table Zeeman 
splitt ing are near their 3c upper limits ijBourke et alJ 
2001). Further more, if GMCs in general are highly flat- 
tened, we ought to observe at least some of them edge- 
on, allowing us to see their sheet-like structure. No such 
sheet-like clouds have been observed. A final problem 
with the sheet-like cloud picture is that a highly flat- 
tened geometry is not consistent with the observation 
that clouds have turbulent energies comparable to their 
gravitational potential energies. In such a cloud, the tur- 
bulence would be strong enough to bring gas out of the 
cloud plane and create a more three-dimensional geome- 
try. 

Another line of observational evidence that clouds are 
magn etically superc r itical comes from statistical indica- 
tors. iPadoan et alJ l)2004|) argue based on simulations 
that the magnetic fields that are at or above equiparti- 
tion with the kinetic energy yield measurably different 
distributions of column density than fields that are be- 
low equipartition. They argue that the observations are 
closer to the sub-equipartition simulations. While there 
is some uncertainty in interpreting simulations of peri- 
odic boxes in the context of real, finite-sized molecular 
clouds, these simulations to provide a strong argument 
for magnetic supercriticality. 

One final problem for magnetically-mediated star for- 
mation theories is that the time required for ambipolar 
diffusion to change a subcritical region into one that is 
supercritical may be consi derably shorter in t urbulent 
media than in static media ijHeitsch et al.ll200^) . Conse- 
quently, the long ambipolar diffusion time invoked to ex- 
plain the low SFR may not apply to GM Cs, which obser- 
vation s indi cate are strongly turbu lent. iLi fc Nakamural 
( 2004) and iNakamura fc Lit l)2005|) perform simulations 
showing that in a two-dimensional geometry, turbulence 
does not enhance ambipolar diffusion enough to make the 
star formation rate too high, but two-dimensional and 
three-dimensional turbulence are very different, so it is 
unclear that their results in this regard are applicable to 
real clouds. Whether magnetic regulation with ambipo- 
lar diffusion is even capable of producing the correct star 
formation time scale in a turbulent medium remains an 
open question. 

7.4. Why Is Q w 1 ? 
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An additional important assumption in our theory is 
that Q ~ 1. While this is well-justified observation- 
allv llQuirklll972t lKennicuttiri989t iMartin fc Kennicuttl 
previous work has also provided a theoretical ex- 
planation, which is part of any complete theory of star 
formation. Theoretically one expects feedback effects to 
prevent Q from straying too far from unity. In ordi- 
nary disk galaxies like the Mil ky Way, su pernovae are 
the likely feedback mechanism (Silk 1997). If Q is too 
low then the star formation rate will increase (equation 
I56|) and the supernova rate will increase as well. This will 
raise the temperature and the velocity dispersion in the 
ISM, increasing Q and reducing the star formation rate. 
If Q becomes too large compared to unity, then gravita- 
tional instability shuts off and molecular clouds cease to 
for m. This is observ ed in the outer parts of disk galax- 
ies (jKennicuttl 1 1 989|) . However, if there is sufficient gas 
present, then without continual supernova stirring the 
gas velocity dispersion will decrease. This will reduce Q, 
causing the star formation rate to rise again. 

In starbursts, the feedback mechanism probably 
changes over from supernovae to ra diation pressure 
(Tho mpson. Quataert. fc Murravl l2005j) . but the effect 
is similar. Low Q values raise the star formation rate, 
which increases the luminosity of the stellar population 
and thereby increases the radiation pressure. This puffs 
up the disk and restores Q«l. If Q is much larger than 
unity, the star formation rate will fall and the disk will 
lose radiation pressure support and begin to collapse, re- 
ducing Q. These mechanisms complete the picture of 
why Q ss 1. 

7.5. Are Molecular Clouds Bound? 

Our analysis also depends on molecular clouds being 
gravitationally bound, virialized structures. If the true 
virial parameter of GMCs in substantially different from 
unity, then SFRg and the overall star formation rate will 
be much greater (for a v j r <C 1) or smaller (for a vir ^> 1) 
than we have estimated. Furthermore, our analysis based 
on the density probability distribution function assumes 
that molecular clouds are gravitationally bound struc- 
tures that live long enough for their density distributions 
to reach statistical equilibrium. If GMCs are largely un- 
bound, or consist of gas that has been co mpressed by 
shock s and that all collapses immediately llElmegreenl 
1 20001 iH artmanm Ballesteros-Paredes 1 fc Berginl 1200 It 
Iciark fc Bonnelll2004 iClark et aUl2005|) . it is not clear 
that the density PDF could reach its equilibrium form be- 
fore the star formation process was complete. We must 
therefore consider whether our assumption of bound, 
virialized clouds is a sound one. 

Observation s indicate that GMC virial parameters are 
close to unity. iMcKee fe Tadj2003[) analyze the CO sur - 
veys of lSolomon et all (|1987^) and iDame et alJ 0986), 
and find that the mean virial parameters for the large 
clouds in their s amples, where mos t stars form, are 1.3 
and 1.4. In M33. IRosolowskv et alJ l)2003|) obtain veloc- 
ity dispersions, masses, and radii for 36 GMCs. From 
their data, we find a mass- weighted mean virial param- 
eter of 1.6. In the nucleus of M64, IRosolowskv fe Blit3 
( 2005aj) find that GMCs are overpressured with respect 
to their environments by at least a factor of 2, indicating 
that they too likely have a V i r ~ 1. Thus, our adopted 
value of a v ir = 1-3 is in good agreement with observa- 



tions, both in the Milky Way and in the disks and nuclei 
of galaxies similar to it. 

That observed virial parameters are all close to unity 
in itself strongly indicates that GMCs are gravitationally 
bound, not held together temporarily by the ram pres- 
sure of turbulent flows in the ISM. There is no reason 
that turbulent flows would create clouds with a v ; r w 1. 
As an example, consider the molecular clumps inside 
GMCs, most of which are created by turbulent flows and 
confined by turbulent pressure rather than self-gravity. 
Most clumps have virial parameters a V ir 3> 1, and they 
have a power-law distribu tion of a va values for a v ; r ^ 1 
ijBertoldi fe McKe"elll992|) . The same is true of molecu- 
lar clouds with masses < 10 4 M l|Hever fe Bruni2004|) . 
While molecules will only form in dense regions of the 
ISM, and for this reason CO surveys are biased towards 
dense gas with low virial parameters, for the UV field of 
our galaxy to be such that we see only clouds that have 
a virial parameters of 1 — 2 requires an unlikely coinci- 
dence. Even if this coincidence could work in the Milky 
Way, it would not explain the observations in M33, where 
the interstellar UV flux could be quite different, and in 
M64, where the density of the gas prevents FUV photons 
from propagating through the ISM at all. 

Another strong argument that suggests GMCs are 
bound is that GMCs have a characteristic mass. In the 
Milky Way, there is a clear upper limit on GMC masses 
of approximately 6 x 10 6 Mq. This limit is not consistent 
with statistically "running out" of clouds at high masses. 
It is a real break in the p ower-law distribution that is ob- 
served at lower masses HMcKee fc Williams! li"997|). The 
mass distributions of GMCs in M33 lEngandola et alJ 
(2003) and M64 (E. Rosolowsky, private communication) 
also exhibit characteristic scales. If GMCs are gravita- 
tionally bound, then the Jeans mass provides a natural 
scale that agrees reasonably well with the observations. 
If GMCs are not bound, however, they cannot have been 
created by gravitational collapse and the Jeans mass is 
therefore irrelevant. Turbulent flows without self-gravity 
are scale-free. If GMCs are unbound they should not 
exhibit any characteristic mass. This prediction of the 
unbound GMC model is inconsistent with the observa- 
tions. One cannot invoke observational selection biases 
to explain this inconsistency, as is done to explain the 
observed values of a v i r . Rendering the Milky Way GMC 
mass distribution consistent with a pure power law would 
require that the Milky Way c ontain ss 100 GMCs with 
masses larger than 6 x 10 6 M Q ijMcKee fc Williamslll997l: 
lMcKeell999l) . There is no plausible way that such a large 
number of very massive clouds could have been missed. 

7.6. Feedback and Cloud Destruction 

Thus far we have omitted any discussion of the ef- 
fects of massive star formation feedback. Obviously mas- 
sive star formation gives rise to HII regions tha t destroy 
molec ular clouds by photoionization and winds. iMatznerl 
( 2002) estimates that this effect limits galactic GMCs to 
converting at most ~ 5 — 10% of their mass into stars 
over their lifetimes. Our justification for neglecting this 
effect hinges on the difference between the star formation 
efficiency, which measures the fraction of gas in a par- 
ticular GMC that is transformed into stars, and the star 
formation rate, which measures the instantaneous rate 
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at which gas is transformed into stars. Feedback from 
massive stars ultimately controls the star formation effi- 
ciency by disrupting a cloud before it can turn most of its 
mass into stars. However, feedback does not change the 
instantaneous star formation rate in the molecular gas 
except indirectly, by driving turbulence in the molecular 
gas and therefore changing the Mach number. Feedback 
only affects the star formation rate by turning molecu- 
lar gas into atomic or ionized gas, thereby reducing the 
amount of molecular gas available to make stars. A thor- 
ough understanding of mechanisms like photoionization 
that regulate the amount of molecular gas available to 
form stars would allow us to calculate /gmc from first 
principles, rather than taking it from observations, and 
would be an important piece of a complete theory of star 
formation. However, our results can stand independently 
of this, since /gmc is directly observable, and our theory 
therefore relies only on direct observables. 

7.7. Turbulent Decay 

The largest single omission from our theory of star 
formation is that it does not address the critical ques- 
tion of what keeps GMCs in virial balance. Simulations 
of both hydrodynamic and magnetohydrodynamic tur- 
bulence in periodic boxes indicate that turbulence de- 
cays on time scales of a si ngle crossing time of the ob- 
ject llMac Low et all Il998t iStone et all 11998 IMac Low! 



ll999HPadoan fc Nordlund 1 1999(1 . The crossing time is 
t cx = 2R/a, where R is the object's radius and a is its 
velocity dispersion. The crossing time and the free-fall 
time are related to the virial parameter by 



5o- 2 i? _ 40 
GM ~ 3^ 



(81) 



where M is the object's mass. In a cloud with our fidu- 
cial value of a V ir = 1-3, t cr = 1.8tg. If the turbulence 
decays substantially in a single crossing time, this means 
that the object should enter free-fall collapse within ~ 2 
free-fall times. If that happened, then a v - n would be- 
come much smaller than unity, and the majority of the 
gas would rapidly turn into stars. That would yield a star 
formation rate far higher than observations allow. Thus, 
GMCs must not be collapsing in this manner. Rapid 
decay of turbulence is also difficult to reconcile with sev- 
eral other observations (see McKee 1999 for a detailed 
discussion). 

There are several possible explanations for the non- 
collapse of GMCs. First is the possibility that turbu- 
lence mB^jioj^_de£a^^s_cjuickly as the simulations indi- 
cate. iCho fc Lazarianl l)2003fl argue that Alfven waves in 
a turbulent magnetized medium cascade from large to 
small scales and decay anisotropically, with modes along 
and perpendicular to the magnetic field having different 
decay rates. Only one mode decays as rapidly as the 
simulations indicate. They argue that the simulations 
performed to date lack t he dynamic range to model this 
effect correctly. Similarlv. lSugimoto. Hanawa. fc Fukudal 
(2004) perform simulations showing that, in a filamen- 
tary cloud geometry, Alfven waves of different polariza- 
tions decay at different rates, with some modes decay- 
ing twice as slowly as earlier simulations indicated. If 
these results from somewhat idealized cases apply to real 
clouds, then GMCs could live for several free-fall times. 



long enough to allow the formation of massive stars that 
could disrupt them rather than letting them collapse en- 
tirely into stars. 

A second possibility is that turbulence in GMCs is 
driven by continual perturbations from outside that are 
strong enough to preve nt the decay of turbulenc e and 
keep GMCs virialized. iKornreich fc Scale! (|2000() sug- 
gest that GMCs will be struck by supernova shock waves 
that maintain cloud turbulence at intervals comparable 
to the free-fall times of large GMCs. However, this source 
of driving is highly stochastic, so it is unclear the shocks 
can truly keep m o st clou ds from collapsing. Furthermore, 
iNakamura et alJ l)2005l) perform numerical studies indi- 
cating that it may not be possible for external shocks 
to drive turbulence in clouds withou t disrupting them 
entirely. iKovama fc Inutsukal l(2002fl suggest that tur- 
bulence is driven by thermal instability at the interface 
between atomic and molecular gas. However, the charac- 
teristic size scale of the disturbances this creates is only 
~ 0.1 pc, so it is unclear that this t urbulence would be 
able t o affect the interiors of GMCs. iPiontek fc Ostrikerl 
(2004) consider thermal plus magnetorotational insta- 
bility in the atomic phase of the ISM, and find that 
magnetic fields allow motions generated at the warm- 
cold interface to drive turbulence far from the interface. 
However, it is unknown if this mechanism would work in 
GMCs. Furthermore, thermal instability offers no clear 
way to explain turbulence in GMCs in galaxies like M64 
where the ISM has no atomic phase, and is not known to 
be thermally bistable as is the atomic ISM in the Milky 
Way. 

A third possible solution to the problem of turbu- 
lent decay is d r iving by fe edback f r om st ar formation. 
iNorman fc Silkl (|1980D and iMcKeei l|1989ft argue based 
on analytic calculations that, for the observed star for- 
mation rate, the rate at which protostellar outflows inject 
energy into their parent clouds is su fficient to balance the 
rate at which turbulence decays. [C)uillen et al. (2005) 
observe protostellar outflow cavities in NGC 1333, and 
estimate that the rate of energy injection from the ob- 
served cavities is sufficient to power t he turbulence of the 
cloud, in agreement with this model. iMatzneJ l|2002f) ar- 
gues that when massive stars are present, turbulent mo- 
tions driven by the overpressure in HII regions are the 
dominant source of energy injection. Matzner estimates 
analytically that the energy injection rate by HII regions 
is sufficient to balance the turbulent decay rate even if 
the decay time is only a crossing time. However, the the- 
ory depends on an efficiency of energy injection by HII 
regions that has only been estimated analytically, and 
ideally should be set by simulations. 

Regardless of the true mechanism, the observations 
show that GMCs cannot be collapsing completely and 
rapidly. The exact mechanism by which the turbulence 
is maintained does not affect our analysis, because, below 
the scale at which it is driven, all turbulence is the same. 
That is why, for example, simulations find a universal 
density probability distribution function independent of 
whether the turbulence is driven or undriven, and re- 
gardless of the random realization of the initial velocity 
field or driving field. Observed GMCs both in the Milky 
Way and in other galaxies are virialized, with turbulence 
balancing gravity, and we have shown here that virial- 
ized, turbulent clouds produce a star formation rate that 
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is consistent with observations. The remaining signifi- 
cant piece of this theory, which we leave for future work, 
is an explanation for how the observed virial balance is 
maintained. 

8. CONCLUSION 

In this work we have attempted to fill in a significant 
missing piece of the overall picture of star formation: 
a quantitative theory that can map the conditions in 
a star forming region into a star formation rate based 
on simple physical principles. Our basic picture is that 
stars form in gravitationally bound, virialized molecular 
clouds. Only 1 — 2% of a cloud is transformed into stars 
in a single free-fall time because in a turbulent virialized 
cloud, most of the gas is in structures that have more 
kinetic energy than gravitational potential energy. Only 
rare, overdense regions are gravitationally bound, and 
the fraction of a cloud's mass in such regions is nearly 
a constant ~ 1% over all virialized clouds. We have for 
the first time computed the collapsing mass fraction di- 
rectly in terms of the Mach number and the virial param- 
eter, the two basic dimensionless numbers that describe 
a star-forming cloud, and shown that the fraction of gas 
in collapsing structures is only a very weak function of 
the Mach number for virialized clouds. The star for- 
mation rate is simply the mass in sufficiently overdense 



structures divided by the cloud free-fall time. Our model 
does not rely on an unknown efficiency of star formation 
or an unknown critical density. The only inputs are the 
physics of turbulence and the virial theorem. 

This prescription correctly predicts the star formation 
rate when we apply it to the observed giant molecular 
clouds in the Milky Way. We also estimate the properties 
of star-forming clouds in other galaxies as a function of 
the rotation speeds and surface densities of various com- 
ponent in those galaxies. We use these estimated cloud 
properties combined with our prediction for the star for- 
mation rate in a cloud to compute galactic-average star 
formation rates, and show that our predictions agree with 
the observed star formation rate in a sample of galaxies 
ranging from normal disks like the Milky Way to star- 
bursts and ULIRGs. Thus, our theory provides a unified 
model capable of explaining the star formation on scales 
from the individual clouds within a galaxy to the entire 
star-forming disk of a starburst or normal disk galaxy. 
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APPENDIX 
ESTIMATING <\> P 

We estimate 4>p by considering cases ra nging from normal disks t o starbursts. In the solar neighborhood, the total 
disk surface density is S tot « 56 M Q pc~ 2 ( Holm berg fc Flvnr]l2004j) . and the gas surface density is S g « 12 M Q pc~ 2 
ijBoulares fc Coxlfl990(l . so / g « 0.21. The total midplane pressure is P ss 3.9 x 10~ 12 dyn cm " 2 , but approximately 
1.9 x 10 -12 dyn cm -2 of this comes from magnetic fields and cosmic rays (jBoulares fc CoxHl990D . Since these permeate 
the molecular clouds and the non-molecular gas equally, they provide no confining pressure on molecular clouds. The 
effective pressure on GMCs in the Milky Way, therefore, is roughly 2 x 10~ 12 dyn cm~ 2 . For the Milky Way solar 
neighborhood values of S to t and S g , we find mp = 0.50. T hus, <t>p ss 2.4 in the s olar n eighborhood. 

At the opposite extreme consider a starburst or ULIRG. iDownes fc Solomon! l)1998l) find that the gas fraction in 
high-surface density starbursts is f g ss 1/3. We cannot directly observe (f> m p in starbursts, but we can estimate it 
based on physical considerations. The reason </> mp < 1 in the Milky Way is that the gas scale height is small compared 
to the stellar scale height. This occurs because the gas comprises a small fractio n of the total s urface density of the 
disk, and because old stars have had a long time to scatter off molecular clouds l|Rafikov 20Q3). In a starburst, the 
gas frac tion is considerably higher and there is no population of old stars that have had a long time to be dynamically 
heated ijDownes fc Solomonlll998T) . We therefore expect that stars and gas will have comparable scale heights, which 
will produce 4> mp ~ 1. This gives <j)p = 3 in starbursts. 

Since (ftp seems roughly constant over a range of environments from the solar neighborhood to extreme starbursts, 
we adopt a constant value of 4>p — 3 throughout our work. The plausible range of variation of <pp is from ~ 1, 
corresponding to a purely gaseous disk, to 4 >p ~ 6, corresponding to a starburst containing only 1/6 gas, the rough 
lower limit in the IDownes fc Solomon! l)1998l) sample. 

Note that, because GMCs occupy a relatively small fraction of the ISM, one might treat them as a pressureless 
component like stars rather than a pressure-contributing component like atomic gas. This would reduce (pp. However, 
since within a GMC the molecular gas does contribute pressure, the product 4>p(j>p must remain unchanged. Thus, if 
one takes a smaller value for <j)p one must use a correspondingly larger value for <j>p. Since our predicted star formation 
rate depends on </>_p</>p, there would be no net change to our predictions. 

ESTIMATING <j>p 

In an environment where the ISM is predominantly atomic, such as the Milky Way, interstellar UV photons dissociate 
H2 and CO that is not sufficient ly shielded. Thus, molecular clouds exist only as the inner parts of atomic-molecular 
complexes l|Elmegreenlll989l Il994|) . Since atomic and molecular hydrogen cannot cool effectively, star formation only 
occurs in the parts of the complexes where CO is present. For Milky Way interstellar UV fluxes, a layer of gas 
wher e C is atomic must provide at least ~ 0.7 mag of extinction to prevent dissociation of CO l|van Dishoeck fc Black 
1988). With such a shielding layer, t he mean pressure in the molecular gas is higher than in the combined atomic and 
molecular complex. iHollimanl l)1995j) estimates (j>p ~ 8, which is consistent with the observed ratios of GMC pressure 
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to ISM pressure in the Milky Way (|BlitzHl99^) . However, there is considerable uncertainty in applying this estimate 
to other galaxies, because it depends on the metallicity of the galactic ISM and the strength of the interstellar UV 
field, both of which vary considerably from galaxy to galaxy. 

For galaxies where the ISM is purely molecular, clouds are not exposed to any external UV flux. In this case, we 
assume that clouds can be described very roughly as polytropic spheres. For a polytropic cloud with P oc r~ kp , 

For an isothermal sphere, kp = 2 so <fi-p — 3. For a cloud with a density profile p oc r -1 , <ftp — 1. We consider these 
extreme limits, and take = 2 as a typical value. This is consistent with observations of GMCs in purely molecular 
galaxies ijRosolowskv fc Blitz 20 0~5aj) . 
We adopt a very rough formula to interpolate between the purely atomic and purely molecular cases: 

0p = 10-8/gmc, (B2) 

where /gmc = S mo i/S g is the molecular gas fraction. One could also have chosen to use a step function approximation 
or simply taken <j?p = 6 as a universal value covering the range from starbursts to ordinary disks. We consider any 
value of 4>p from 2 to 10 reasonable, although a value of 2 is implausible for a galaxy with a great deal of atomic gas, 
and a value of 10 is implausible for a galaxy that is entirely molecular. 
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